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Abstract 

A dynamic procedure for the Lagrangian Averaged Navier-Stokes-a (LANS-a) equations is devel- 
oped where the variation in the parameter a in the direction of anisotropy is determined in a 
self-consistent way from data contained in the simulation itself. In order to derive this model, the 
incompressible Navier-Stokes equations are Helmholtz-filtered at the grid and a test filter levels. A 
Germano type identity is derived by comparing the filtered subgrid scale stress terms with those 
given in the LANS-a equations. Assuming constant a in homogenous directions of the flow and 
averaging in these directions, results in a nonlinear equation for the parameter a, which determines 
the variation of a in the non-homogeneous directions or in time. Consequently, the parameter a 
is calculated during the simulation instead of a pre-defined value. The dynamic model is initially 
tested in forced and decaying isotropic turbulent flows where a is constant in space but it is al- 
lowed to vary in time. It is observed that by using the dynamic LANS-a procedure a more accurate 
simulation of the isotropic homogeneous turbulence is achieved. The energy spectra and the total 
kinetic energy decay are captured more accurately as compared with the LANS-a simulations us- 
ing a fixed a. In order to evaluate the applicability of the dynamic LANS-a model in anisotropic 
turbulence, a priori test of a turbulent channel fiow is performed. It is found that the parameter a 
changes in the wall normal direction. Near a solid wall, the length scale a is seen to depend on the 
distance from the wall with a vanishing value at the wall. On the other hand, away from the wall, 
where the turbulence is more isotropic, a approaches an almost constant value. Furthermore, the 
behavior of the subgrid scale stresses in the near wall region is captured accurately by the dynamic 
LANS-a model. The dynamic LANS-a model has the potential to extend the applicability of the 
LANS-a equations to more complicated anisotropic fiows. 
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I Introduction 



Turbulent flows play a an important role in many areas of engineering fluid mechanics as well as 
atmospheric and oceanic flows. Accurate simulation of a turbulent flow requires that the energetics 
of the large scale energy containing eddies, dissipative small scales, and inter-scale interactions to be 
accounted for. In direct numerical simulations (DNS) all the involved scales are directly calculated. 
DNS is believed to provide the most comprehensive representation of the governing equations of 
fluid flows; the so-called Navier-Stokes (NS) equations. Owing to the very high Reynolds numbers 
encountered in most problems of interest, the disparity between the large scales and small scales, 
which represents the computational size of the problem, rapidly grows with the Reynolds number. 
Consequently, DNS can resolve only a small fraction of the turbulent activity for high Reynolds 
number flows. 

While the direct numerical simulation of most engineering flows seems unlikely in near future, 
turbulence modeling could provide qualitative and in some cases quantitative measures for many 
applications. Large Eddy Simulations (LES) and the Reynolds Averaged Navier-Stokes Equations 
(RANS) are among the numerical techniques to reduce the computational intensity of turbulent 
calculations. In LES, the dynamics of the large turbulence length scales are simulated accurately 
and the small scales are modeled. The vast majority of contemporary LES make use of eddy- 
viscosity based Subgrid-Scale (SGS) models in conjunction with the spatially-averaged (filtered) 
Navier-Stokes Equations. In this approach, the effect of the unresolved turbulence is modeled as 
an effective increase in the molecular viscosity. On the other hand, RANS models are obtained by 
time averaging the Navier-Stokes equations. In this case most of the unsteadiness is averaged out. 
Consequently, the time mean quantities are calculated while the faster scale dynamics are modeled. 
RANS simulations are often more affordable than LES, however, their accuracy is somewhat limited 
in many applications [24]. 

More recently, Holm, Marsden and their coworkers [12] introduced a Lagrangian averaging 
technique for the mean motion of ideal incompressible flows. Figure 1 contrasts the derivation 
of LES, RANS, and the Lagrangian Averaged Navier-Stokes-a (LANS-a) equations. Unlike the 
traditional averaging or filtering approach used for both RANS and LES, where the Navier-Stokes 
equations are averaged or spatially filtered, the Lagrangian averaging approach is based on averaging 
at the level of the variational principle. In the isotropic Lagrangian Averaged Euler-a (LAE-a) 
equations, fluctuations smaller than a specified scale a are averaged at the level of the flow maps 
[3] . Mean fluid dynamics are derived by applying an averaging procedure to the action principle of 
the Euler equations. As shown in Figure 1, both the Euler and the Navier-Stokes equations can be 
derived in this manner (see Marsden & Ratiu [18] for a variational derivation of the Euler equations). 
The usual Reynolds Averaged Navier-Stokes (RANS) or LES equations are then obtained through 
the subsequent application of either a temporal or spatial average. The critical difference with the 
Lagrangian averaging procedure is that the Lagrangian (kinetic energy minus potential energy) is 
averaged "prior to the application of Hamilton principle and a closure assumption is applied at this 
stage. This procedure results in either the Lagrangian averaged Euler Equations (LAE-a) ^ or the 
Lagrangian averaged Navier-Stokes Equations (LANS-a), depending on whether or not a random 
walk component is added in order to produce a true molecular diffusion term. Since the Hamilton 
principle is applied after the Lagrangian averaging is performed, all the geometrical properties (e.g. 
invariants) of the inviscid dynamics are retained even in the presence of the model terms which arise 

^In this nomenclature, a is used to denote the filtering scale (i.e. the simulation faithfully represents motions on 
a scale larger than q). 
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Figure 1: Derivation of the averaged flow equations. 



from the closure assumption [12; 11; 10]. For instance, LAE equations posses a Kelvin circulation 
theorem. Thus it is potentially possible to model the transfer of energy to the unresolved scales 
without an incorrect attenuation of quantities such as resolved circulation. This is an important 
distinction for many engineering and geophysical flows where the accurate prediction of circulation 
is highly desirable. 

Numerical simulations by Chen et al [5] and Mohseni et al [22] showed the capability of the 
LANS-a equations in simulating isotropic homogenous turbulence. However, most engineering and 
geophysical flows of interest are often anisotropic. For example, due to rapid damping of turbulent 
fluctuations in the vicinity of a wall, the application of the isotropic LANS-a equations with a 
constant a is not appropriate for long term calculations. In order to capture the correct behavior 
in such systems the parameter a must be spatially or /and temporally varied in the direction of 
anisotropy [4], i.e. wall normal direction. There has been some attempt (with limited success) in 
order to remedy this problem. A successful dynamic LANS-a model is yet to be formulated and 
tested. There are at least two approaches to anisotropy in the LANS-a equations: 

(i) To derive a set of anisotropic LANS-a equations. See alternative derivations in [10; 19]. 

(ii) Use the isotropic LANS-a equations, but with a variable a to compensate for the anisotropy. 

At this point much more work must be done on the anisotropic LANS-a equations before they can 
be applied to practical problems. The second approach listed above is what will be explored in this 
study. 

This paper is organized as follows: The isotropic LANS-a equations and some of their main 
features are summarized in section 2. A dynamic LANS-a approach is proposed in section 3 where 
the variation in the parameter a in the direction of anisotropy is determined in a self-consistent way 
from the data contained in the simulation itself. Our approach will be developed in the same spirit as 
the dynamic modeling procedure for conventional LES [8; 9; 20; 13] which has achieved widespread 
use as very effective means of estimating model parameters as a function of space and time as the 
simulation progresses. The incompressible Navier-Stokes equations are Helmholtz-filtered at the 
grid and a test filter levels. A Germano type identity is derived by comparing the filtered subgrid 
scale stress terms with those given in the LANS-a equations. Considering a constant value of a and 
averaging in the homogenous directions of the flow results in a nonlinear equation for the parameter 
a, which determines the variation of a in the non-homogeneous directions or time. This nonlinear 
equation is solved by an iterative technique. Consequently, the parameter a is calculated during 
the simulation instead of a fixed and pre-defined value. 

Numerical techniques for simulating the dynamic LANS-a model in this study are described 
in section 4. The performance of the dynamic LANS-a model in simulating forced and decaying 
isotropic homogeneous turbulent flows are considered in section 5. In statistically equilibrated 
forced turbulence, the parameter a should remain constant in time and space. In decaying isotropic 
turbulence, the parameter a could change in time as the integral scales of the turbulent flow changes. 
In order to demonstrate the applicability of the dynamic LANS-a model of this study in anisotropic 
flows, a priori test of turbulent channel flows are also performed in section 5. Concluding results 
are presented in section 6. 



2 The Isotropic LANS-a Equations 
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The incompressible isotropic LANS-ct equations for the large scale velocity u are given by (see [12] 
for a derivation) 

du 1 

— + (n-V)u = -Vp+— Am + V-t(u), (1) 
at Re 

V-u = 0, (2) 

where t{u) is the subgrid stress tensor defined as [3] 

r(n) = -0^(1 - a^A)"^ [Vn • Vn^ - Vu^ • Vn + Vu • Vn + Vu^ • Vn'^] . (3) 

The subgrid scale stress t(u) is in fact the momentum flux of the large scales caused by the action 
of smaller, unresolved scales. Here a is a constant length scale introduced during the averaging 
process. Note that for vanishing parameter a the NS equations will be recovered. 
The LANS-a equations can be represented equivalently by 



dv 1 

— — \- (u ■ V)v + VjVuj = —VP -\ At;, where Vi is defined as v = u — a'^Au. (4) 

at J J 



The modified pressure P in these equations is determined, as usual, from the incompressibility 
condition: V • u = and V • t; = 0. 

One interpretation for the equations (1) is that they are obtained by averaging the Euler equa- 
tions in Lagrangian representation over rapid fluctuations whose scale are of order a. In this 
respect, one can show that the Lagrangian averaged Euler equations can be regarded as geodesic 
equations for the metric on the volume preserving diffeomorphism group, as Arnold [2] did with 
the L2 metric for the Euler equations. Note that in calculating the SGS stress t{u) in equation 
(3) one needs to calculate the inverse of the Helmholtz operator (1 — a^A), which implies the need 
to solve a Poisson equation. While efficient numerical treatment of the Poisson equation, or its 
possible elimination through rational approximation will be a focus of a future publication, we note, 
in passing, that the inverse of the Helmholtz operator can be expanded in a to higher orders of the 
Laplacian operator as shown in below 

(1 - a^A)-! = 1 + a^A + a^A^ + • • • . 

Consequently, solving a Poission equation for inverting the Helmholtz operator could be avoided. 

It is interesting to note that the Lagrangian averaging technique preserves the Hamiltonian 
structure of the governing equations in the inviscid limit while the effects of small scales on the 
macroscopic features of large scale are taken into account in a conservative manner. The Hamil- 
tonian and Lagrangian formulations of ideal fluids are both basic and useful. These formulations 
are part of a more general framework of geometric mechanics, which plays a vital role in the de- 
velopment of new continuum models suited for computation, as well as numerical algorithms that 
preserve structure at the discrete level. In recent years the geometric approach to fluid mechanics 
has been quite successful. Geometrical methods provide a framework for the study of nonlinear 
stability [1], variational integrators [14; 16], statistical equilibrium theory [6; 21], and many other 
interesting topics in fluid dynamics. The Lagrangian averaged Navier-Stokes-a uses ideas from 
geometric mechanics and offers a theoretically and computationally attractive approach to the 
turbulence closure problem. 



S Derivation of a Dynamic LANS-a Model 



The LANS-a equations for the large scale velocity u are given by equations (1), where t{u) is 
the subgrid stress tensor defined in (3). This set of equations for a is similar to the grid filtered 
equation in the dynamic LES. In analogy with the dynamic LES one can obtain an equation for 
the filtering length scale, a, by filtering the Navier-Stokes equations 

duj ^ ^ duj _ dp ^ I d'^Uj 
dt dxj dxi Re dxjdxj ' 

with the Helmholtz related filters 



u 
u 



to obtain 



(l-a^A)-^^, grid filter, (6) 

(1 - a2A)~i(l - a^A)"^^, test filter, (7) 



duj ^ dujUj _ dp ^ I d'^Uj dnj 
dt dxj dxi Re dxjdxj dxj ' 

dui ^ dujUj _ _dp_ _^ 1 d'^Uj _ dTjj 



dt dxj dxi Re dxjdxj dxj 



where 



J- ij — Uj^Uj ttittj. 

Using an idea similar to Germano identity [8] , we define 

Lij = Tij — fij = UiUj — UiUj, (10) 

where the subgrid scale stresses under two filtering actions can be modeled by the LANS-a subgrid 
term in equation (3). Therefore, 

Tij = a\l-a'^A)-^Mij, (11) 



where 



Ti^=a\l-a^A)-^Nij, (12) 



dxk dxk dxi dxj dx^ dxj dx^ dxi ' 
dui duj dukduk dui duk du^ duk 

Nij = TT^-T^T^ V- h 



dxk dxk dxi dxj dxk dxj dxk dxi 
Combining equations (10)-(12), one obtains 

Lij = p^a\l - P^a'^A)-^Nij - a\l - P^a^A)~\l - a^A^^Mij, (13) 

or 

Lij = a^iP^Ni, - Mi,), (14) 



1/a 1/a 
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Figure 2: The positions of grid and test filter scales on the turbulent kinetic energy spectrum. 

where j3 = a/a. Multiplying both sides of the above equation by 5ij, to yield 

UjSij = a\(3^Nij - Mij)Sij. (15) 
Taking spatial averaging of both sides of the above equation in homogenous directions, one obtains 

a2 = (LijSij) ^^g^ 

where 

^ 1 / dui duj 
^■^ 2 \dxj dxi 

The denominator in equation (16) could approach zero, where it creates a singularity. In dynamic 
LES, Lilly [17] used a least square approach to eliminate the singularity in Germano's model. By 
a similar least square approach a nonlinear equation for a could be found as 

«2 ^ ^ (L.,(/32iV.,-4,)) ^^^^ 

((/32iV,,-M,j)(/3'%- Ai))' 

which does not have the singularity problem as in equation (16). This is a nonlinear equation for 
a. All the quantities in equation (17) can be calculated during a LANS-a simulation. Therefore, 
equation (17) provides a nonlinear equation for dynamically calculating the value of a during the 
simulation. 

At this point the potential values for the free parameter (3 are required. Writing the grid and 
test filters in equations (6) and (7) in the Fourier space, one obtains 

and 

{l + P^a^k^){l + a^k^)^ l + iP^ + l)a^k^ 1 + a^k^ k ^ oo, ^yj 

where (•) stands for variables in the Fourier space, k is the wavenumber, and a corresponds to 
filter scale for the test filter. Since a = ylT/J^a > a, one can realize that as long as /3 > 0, 



&ie test filter have a larger filter scale than the grid filter. Figure 2 shows the relative positions of 
the grid filter scale a and the test filter scale a on a schematic of the energy spectrum for a high 
Reynolds number flow. In order to accurately model the subgrid scale stress, both the grid filter 
and the test filter scales must be located in the inertial sub-range of the energy spectrum. It should 
be pointed out that the iterative calculation required in equation (17) does not require new flow 
field calculations, and the iteration at each time step is carried out using the existing flow field at 
the same time step. Similar to the dynamic LES model, the present dynamic LANS-a model has 
a free parameter /3, which is related to the characteristic length scale of the grid and test filters. 

The dynamic a model given in equation (17) is designed to capture the length scale variations 
in space and time. Aside from the isotropic homogenous turbulent flows, it is well suited for 
anisotropic flows such as wall bounded turbulence or mixing flow turbulence, where the turbulence 
length scales could change in space or in time. In cases where there are directions of homogeneity, 
such as the streamwise and spanwise direction in a channel flow, one can average the parameter 
a over the homogeneous directions. In a more general situation, we expect to replace the plane 
average, used in the channel flow, by an appropriate local spatial or time averaging scheme. For 
isotropic homogenous turbulence, a is regarded as a constant in space and changes only in time. 

4 Numerical Method 

The dynamic procedure in this study is initially tested for forced and decaying isotropic turbulence 
where the parameter a is constant over the computational domain, but can vary in time. Further- 
more, a "priori test of the dynamic LANS-a procedure in a turbulent channel flow is investigated. 
In this section the numerical technique for solving the governing equations are summarized. 

Isotropic homogeneous turbulence. The computations are performed in a periodic cubic 
box of side 27r. A standard parallel pesudospectral scheme with periodic boundary conditions 
are employed. The spatial derivatives are calculated in the Fourier domain, while the nolinear 
convective terms are computed in the physical space. A fourth order Runge-Kutta scheme is 
implemented to advance the flow field in time. The two third rule is used in order to eliminate the 
aliasing errors. Therefore, the upper one third of the wave modes are discarded at each stage of the 
Runge-Kutta scheme. The initial velocity field for each case was divergence free and constructed 
to generate an energy spectrum of the form 

E{k) ~ k^exv\-2(klkpf\. 

The value of hp corresponds to the peak in the energy spectrum. The initial pressure fluctuations 
were obtained from the solution of a Possion equation. 

Turbulent flow in a channel. DNS data from del Alamo and J. Jimenez [7] are employed for 
the a priori test. The computational domain in this case, normalized based on the half channel 
height, is spanned Svr in the streamwise and 47r in the spanwise directions. The spatial derivatives 
are calculated by the pesudospectral method in streamwise and spanwise directions and by the 
Chebychev-tau technique in the wall normal direction. Similar computational techniques have 
successfully been used for the DNS of channel flows by Kim et al [15] and Moser et al [23]. Grid 
and test filters of Helmholtz types are applied in both streamwise and spanwise directions, while 
no explicit filters are applied in the wall normal direction, a is assumed to be constant in the 



homogenous directions, i.e. the streamwise and spanwise directions, in order to solve the nonhnea? 
equations (16) or (17). These equations are solved by an iterative technique. Since both the 
mean flow and the flow perturbations vanish at the wall, singular behavior might occur in these 
equations. This can be easily fixed by starting the a priori test a few grid points away from the 
wall. In actual simulation of the dynamic LANS-a equations, one can explicitly put a to zero below 
in the immediate vicinity of a wall when the value of a drops below a threshold. The converged a 
values at each point is used as an initial value for the iteration at the next grid layer. 



5 Results and Discussions 



Capabilities of the dynamic LANS-a model of the previous sections are examined in both isotropic 
and anisotropic turbulent flows. In isotropic homogeneous turbulence the parameter a is constant 
in space but allowed to vary in time. Results of the dynamic model is compared with the isotropic 
LANS-a simulations with a constant a and with the DNS data. 



Decaying isotropic homogenous turbulence simulations. DNS of a decaying isotropic ho- 
mogenous turbulence with initial Taylor Reynolds number of Rex = 72 (corresponding to a com- 
putational Reynolds number Re = 300) is performed to be used as a test case. The initial energy 
spectrum is peaked at kp = 4. The isotropic LANS-a and the dynamic LANS-a simulations are 
calculated for both 64^ (corresponds to 48'^ after dealiasing) and 48'^ (corresponds to 32^ after 
dealiasing) resolutions, and direct numerical simulations are performed for 128^ (corresponds to 
85^ after dealiasing). The eddy turn over time for this case is found to be r = 0.9. Figure 3 shows 
the time evolution of a for (3 = 0.8,0.9, 1, and 1.2. The values of a experience a sharp decrease 
from its initial value during the first eddy turn over time. However, it quickly settle down toward 
a much slower varying value. Slight changes in a value after the first eddy turn over time could be 
traced back to flattening of the energy spectrum as the turbulence decays. 
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Figure 3: Evolution of a for different f3 for a decaying isotropic turbulence at Rex = 72 and r = 0.9. 
Grid resolution (a) 48^ (b) 32^. 



The energy spectra at two different times are shown in Figure 4, and the total kinetic energy 
decay are shown in Figure 5. While a slight dependency on the value of /3 is observed, in general, the 
energy spectrum at various times and the total kinetic energy decay are captured nicely. Mohseni 
et al [22] demonstrated that in order to accurately simulate a turbulent flow with the LANS-a 
equations, the value of a should be somewhere, perhaps one decade lower than the peak of the 
energy spectra toward the grid resolution. Careful considerations of Figures 3 and 4 reveal that the 
dynamic LANS-a model of this study satisfies this criteria for all (3 values. In general, one expects 
that the value of a to be in the inertial range of the energy spectra in order to correctly capture 
the dynamics of the large scales. As illustrated in Figures 4 and 5, it is evident that the dynamic 
LANS-a model provides a better estimate of the total kinetic energy decay and the energy spectra 
over similar simulations with fixed a calculations. 
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Figure 4: Energy spectra of the DNS, dynamic LANS-a, and LANS-a with fixed a simulations of 
a decaying isotropic turbulence at Re\ = 72 and r = 0.9. Grid resolution (a) 48^, (b) 32"^. 



Forced isotropic homogenous turbulence simulations. Forced isotropic turbulence is one 
of the most idealized and extensively simulated turbulent flows. The numerical forcing of a tur- 
bulent flow is usually referred to the artificial addition of energy at the large scales in a numerical 
simulation. Statistical equilibrium is signified by the balance between the input of kinetic energy 
through the forcing and its output through the viscous dissipation. In this study, we adopted a 
forcing method used in Chen et al [5] and Mosheni et al [22] where the wave modes in a spherical 
shell |i^| = /cq of certain width are forced in such a way that the forcing spectrum follows the 
Kolmogorov —5/3 scaling law, that is 



= ^^^A;-5/3. (20) 



Here fi and Ui are Fourier transforms of the forcing vector /j and velocity Ui, N is the number of 
forced wave modes, and 6 controls the injection rate of energy at the large scales. This particular 




Figure 5: Total kinetic energy decay of the DNS, dynamic LANS-a, and LANS-a with fixed a 
simulations of a decaying isotropic turbulence at Re.\ = 72 and r = 0.9. Grid resolution (a) 48^, 
(b) 323. 



forcing technique enforces the energy cascade in the inertial range starting from the first wave 
mode. In this simulations we choose /cq = 2 and = 0.1. The initial Taylor Reynolds number is 
Re\ = 415 and the initial energy spectrum is peaked at = 1, while the eddy turn over time is 
found to be T = 3.8. The grid resolution for simulations using the dynamic LANS-a equations and 
the LANS-a equations with fixed a is 64^, while the DNS data is performed at a grid resolution of 
128^ before deahasing. 

Figure 6 shows the evolution of a for f3 = 0.8 and 1.0. Similar to the decaying turbulence, a 
sharp decrease in the value of a is observed over the first eddy turn over time, where the values of 
a settles down toward a constant value. This corresponds to an statistically equilibrated state. As 
expected, the final value of a is in the inertial range of the energy spectrum. 

Figure 7 shows the energy spectrum at t = 5.8r for f3 = 0.8 and 1.0. An inertial subrange 
with ~ fc~^/3 energy spectrum is evident in the dynamic LANS-a simulations. The results of 
the dynamic LANS-a simulations are compared with the DNS and the LANS-a simulations with 
a = 0.2. The energy spectra of the dynamic LANS-a simulations for /3 = 0.8 and 1.0 show a better 
agreement with the DNS data than the energy spectra for a LANS-a simulation with a constant a. 

A priori test of a turbulent channel flow. A priori test of the dynamic LANS-a model is 
carried out in order to determine the accuracy of the model in predicting the SGS stresses and 
the energy dissipation rates in a wall bounded flow. The tests are performed on a DNS data of 
del Alamo and Jimenez [7] for a turbulent channel flow. The turbulence Reynolds number, based 
on the wall friction velocity, is Rcr = 550 and the computational grid is 1536 x 257 x 1536 in 
the streamwise, wall normal, and spanwise directions, respectively. After dealiasing the physically 
relevant part of the computational domain reduces to 1024 x 257 x 1023. The mean velocity proflle, 
non-dimensionalized by the wall-shear velocity, is depicted in Figure 8(a), where a log layer from 
« 80 to 220 is observed. Figure 8(b) shows the turbulence intensity proflles from the wall to the 
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Figure 6: Evolution of a for different /? in the forced turbulence case with Ke\ = 415 and r = 3.8. 

middle of the channel in global coordinate which is normalized by half channel height 6. Maximum 
turbulence intensities in all directions are located in the wall layer. 

Figure 9 shows the variation of a with the distance from the wall in both global and wall 
coordinates for (3 = 0.8, 1.0, and 1.2. As demonstrated in Figure 9(b), a values experience a sharp 
increase in the vicinity of the wall up to = 100. This region of sharp increase in the value of a 
contains both the viscous sublayer and the buffer layer. Diminishing values of a is observed as one 
approaches the wall. This is consistent with theoretical expectations that the NS equations ought 
to be recovered in the laminar layer at the wall. The profile of a in the vicinity of the wall shows 
minimal dependency on (5. Away from the wall and beyond = 100, a shows little variation 
across the channel. One can argue that the dynamic LANS-q equations in this case divides the 
flow into two distinct regions: a near wall region that includes both the viscous sublayer and the 
buffer layer where a is a function of the distance from the wall, and a constant a region which 
includes the log layer and the outer layer. In the near wall region a keeps an almost log relation 
with the distance from the wall in wall units. In summary, one can argue that in wall bounded 
flows, the isotropic LANS-a calculations could be used with a constant a beyond y~^ = 100 and 
with a logarithmic relation in the near region. This projection requires further investigation in 
LANS-a calculations. 

Similar to the dynamic LES, one expects that the accuracy of the dynamic LANS-a model to 
depend on its capability of accurately modeling the subgrid scale stresses. The modeled and the 
exact SGS stresses in this flow are shown in Figure 10 for the isotropic component < rn > and 
Figure 11 for the shear stress component < ri2 >, where < • > stands for averaging in streamwise 
and spanwise directions. The general trend of the SGS stresses are captured in the dynamic LANS- 
a model without any ad hoc damping function. Good agreement between the modeled and the 
exact SGS stresses in the near wall region are observed. The SGS stresses vanish at the wall and 
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Figure 7: Energy spectra at t = 5.8r for the DNS, dynamic LANS-a and LANS-a with fixed a 
simulations of a forced isotropic turbulence with Rex = 415 and r = 3.8. 

in the middle of the channel with a maximum value within the wall layer. The exact and modeled 
dissipation < e > are compared in Figure 12. Both the SGS stresses and the modeled dissipations 
are effectively insensitive to variation in /3. 

6 Conclusions 

A dynamic LANS-a model is proposed where the variation in the parameter a in the direction 
of anisotropy is determined in a self-consistent way from data contained in the simulation itself. 
The model results in a nonlinear equation for a. Numerical experiments for decaying and forced 
homogenous isotropic turbulence are performed using the dynamic LANS-a model. The simulation 
results in both cases show an improvement over the LANS-a simulations with a fixed a. 

A priori test of the dynamic LANS-a model in a channel flow is carried out, where good 
agreement between the dynamic LANS-a predictions and the DNS data is observed. The parameter 
a is found to rapidly change in the wall normal direction in the vicinity of the wall. Near the solid 
wall, the length scale a shows a logarithmic dependence on the wall normal direction in wall units. 
Away from the wall, and in the middle of the channel, a approaches an essentially constant value. 
As a result, the turbulent flow is divided into two regions: a constant a region away from the 
wall and a near wall region. In the near wall region, a keeps an almost logarithmic relation with 
the distance from the wall. Consequently, one can argue that in wall bounded flows, the isotropic 
LANS-a calculations could be used with a constant a beyond y'^ = 100 and with a logarithmic 
relation in the near region. These results indicates a promising application of the dynamic LANS-a 
model in wall bounded turbulent flow simulations. 
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Figure 11: The averaged subgrid scale shear stress < ri2 > in global units. 
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Figure 12: The averaged dissipation in global units. 



